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Highly energetic, variable and distant sources such as Active Galactic Nuclei provide a good 
opportunity to evaluate effects due to the emission and the propagation of high energy photons. 
^q' In this note, a study of possible energy-dependent time-lags with PKS 2155-304 light curve as 

measured by H.E.S.S. in July 2006 is presented. These time-lags could either come from the 
emission processes or also sign a Lorentz Symmetry breaking as predicted in some Quantum 
Gravity models. A Cross-Correlation function and a Wavelet Transform were used to measure 
the time-lags. The 95% Confidence Limit on the Quantum Gravity energy scale based on the 
statistical and systematic error evaluation was found to be 7 x 10 17 GeV considering a linear 
correction in the standard photon dispersion relations and assuming that emission-induced 
£^ time-lags are negligible. For now, this limit is the best ever obtained with a blazar. 
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1 Introduction 

Quantum Gravity phenomenology^ has known a growing interest in the past decade, especially 
since the time it was argued that the quantum nature of space-time could have measurable effects 
on photon propagation over large distances. Namely, the quantum nature of the space-time at 
the Planck scale could induce a dependance of the light group velocity with the energy of the 
photons'^! 

It is generally assumed that the light velocity is modified following to 
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at the second order, where Ep = 1.22 x 10 19 GeV is the Planck energy and where £ and £ are 
free parameters which need to be determined. In the following, only the linear correction will 
be considered and a limit on -Bqg = |£| _1 Ep will be set. 

As Ep is very large, the modification of the velocity of the photons is expected to be tiny. 
However, as it is related to the nature of the space-time, it was proposed^ that these tiny 
effects could add up over long distances and lead to detectable time-lags between photons of 
different energies, assuming these photons were emitted at the same time. Taking into account 
the expansion of the Universe^, the time-lag At is obtained from the dispersion relations (pQ): 
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Figure 1: Light curve of the PKS2155-304 flare during the night of July 28 in 2006. Top: 200-800 GeV. Bottom: 
> 800 GeV. The original data points (in black) are binned in two-minute time intervals. The zero time point is 
set to MJD 53944.02. Gray points show the oversampled light curve, for which the two-minute bins are shifted 

in units of five seconds. 



in the case of the linear correction and where f2 m , and Ho are parameters of the Cosmological 
Standard Model (Q m = 0.3, A = 0.7 and H = 70 km s" 1 Mpc" 1 ). AE is the energy difference 
of the photons. The time-lag may decrease or increase with AE depending on the model in 
consideration. The integral term increases with the redshift and takes into account the fact that 
two photons travelling with two different speeds do not take the same time to cross the Universe 
and then do not see the same expansion. 

The goal of the 'time of flight' studies is to measure the time-lags, which should be energy 
dependant. For this, two kind of variable extragalactic sources may be considered: the Gamma- 
Ray Bursts (GRBs) and the Active Galactic Nuclei (AGNs). However, both these sources could 
introduce intrinsic time lags in the measurements. Therefore, it would be necessary to study the 
possible effects as a function of the redshift of the source. When this is not possible, the intrinsic 
(or source) effects are often assumed to be negligible. In the present study, this hypothesis was 
adopted. 



2 Data and methods 

In this note (more details are available elsewhere ®) , the search for propagation effects with 
the blazar PKS 2155-304 is described. During the night of the 28th of July 2006, the H.E.S.S. 
experiment detected^ an exceptional flare of this source, located at z = 0.116, with a high flux 
(10 000 photons recorded in 1.5 hours) and a high variability (rise and fall times of ~200 s). 
The over-sampled light curve of the flare is shown in Fig. [1] in two different energy bands. 

To measure the time lags between photons in two different light curves, two independant 
analyses were carried out, using two different methods: 

• The position of the maximum of the Modified Cross Correlation Function^ (MCCF) gives 
directly the value of the time lag. This method was applied to the oversampled light curves 
of Fig. [T] in the energy bands 200-800 GeV and > 800 GeV (Fig. [2j left) which corresponds 
to < AE > ~1 TeV; 

• Following Ellis et al. the Continuous Wavelet Transform^ (CWT) was used to locate 
the extrema of the light curves with great precision. An extremum of the low energy 
band was associated with an extremum in the high energy band to form a pair. Selection 
criteria were applied to reject fake extrema. The energy bands 210-250 GeV and > 600 
GeV were used, with a bin width of one minute and no oversampling of the light curves. 
These energy bands give a mean AE of ~0.92 TeV. 




Figure 2: Left: the MCCF obtained with the two light curves of Fig. \T\ Right: the MCCF peak distribution 
(CCPD) obtained with 10000 realizations of the light curves. 
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Figure 3: Mean of the MCCF peak distribution 
(CCPD) as a function of the injected dispersion. 
Each CCPD is obtained from ten thousand simu- 
lated light curves. The points have been shifted by 
the mean value of the CCPD of the original data 
(Fig. [2} • The dashed line shows the linear response 
function. 
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Figure 4: The dispersion measured as a function 
of the redshift. The results of the Whipple and the 
MAGIC experiments are compared with the present 
H.E.S.S. result. 



As shown in Fig. [5] (left), the MCCF was fitted with the sum of a gaussian and a first degree 
polynomial to determine the position of the maximum. It gave r pea k = 20 s. In order to evaluate 
the uncertainties on this result, 10000 light curves were simulated in each energy band varying 
the flux within the error bars. For each pair of light curves, the MCCF was computed and its 
maximum was filled into the distribution shown in Fig. [2] (right). This distribution is slightly 
asymmetric, with a mean of 25 s and an RMS of 28 s. Another test, performed by injecting 
a dispersion in the data (Fig. [3]), showed no significant deviation of the measured lag. Then, 
as the time-lag obtained was compatible with zero, a 95% confidence upper limit on the linear 
dispersion of 73 s/TeV was found. 

With the CWT method, two pairs of extrema were obtained giving a mean time delay of 27 
seconds. A method similar to the one used for the CCF was applied to determine the errors, 
which were found to be in a range between 30 and 36 seconds. A 95% confidence limit of 
100 s/TeV was obtained for the linear correction. 

3 Results and discussion 

The results obtained with the two methods are summarized in Table [T] (see next page) . The 
MCCF leads to a limit of£ QG > 7.2 xlO 17 GeV. The CWT gives a lower limit of E QG > 5.2xl0 17 
GeV, mainly due to a larger measured time-lag and a lower value of < AE >. 



Table 1: The results obtained with the two methods in the case of a linear dispersion in energy. 



Method 


Energy bands 


<AE> 


(Ai/A£) 95% CL 


-^QG 95% CL 


MCCF 


200 GeV < E < 800 GeV 
and E > 800 GeV 


1.02 TeV 


< 73 s/TeV 


> 7.2 x 10 17 GeV 


CWT 


210 GeV < E < 250 GeV 
and E > 600 GeV 


0.92 TeV 


< 100 s/TeV 


> 5.2 x 10 17 GeV 



Fig. H] (previous page) shows all the results obtained so far with AGNs. The Whipple 
collaboration set a limit of 4 x 10 16 GeV using a flare of Mkn 421 (z = 0.031) in 199612]. More 
recently, the MAGIC collaboration obtained a limit of 2.6 x 10 17 GeV with a flare of Mkn 501 (z 
= 0.034)12. The results obtained with H.E.S.S. are more constraining due to the fact that (i) 
PKS 2155-304 is almost four times more distant than Mkn 421 and Mkn 501, (ii) the statistics 
are higher by a factor of ten. 

As mentionned above, GRBs are other good candidates for time of flight studies. Population 
studies have already been carried out^^l and lead to limits of the order of 10 16 GeV. The best 
limit so far is Eqq > 2 x 10 18 GeV. It has been obtained by the Fermi experiment 1^1 with 
GRB 080916C (z = 4.35) using both GBM and LAT data. However, the time lag was obtained 
comparing the arrival time of the highest energy photon (13.2 GeV) with the trigger time, 
leading to At = 16.5 s. A more detailed analysis including error calibration would be necessary. 

As far as the search for the Quantum Gravity propagation effect is concerned, the GRBs 
and the AGNs are complementary. AGN flares can be detected with high statistics with ground- 
based gamma-ray telescopes, which give large values of AE ~ 1 TeV. However, the absorption 
of the high energy photons by the extra-galactic background light limits the distance of the 
observed objects. On the other hand, GRBs are easilly detected by satellite experiments at 
very high redshifts (up to ~6) and up to a few hundreds of GeV (AE ~ 10 GeV). Till now, no 
significant result has been obtained with either AGNs or GRBs for the linear and quadratic terms 
in the photon dispersion relations. Further observations of both a high number of AGN flares 
and GRBs will be necessary to give robust conclusions on possible propagation effects. Present 
and future experiments such as Fermi, CTA or 

AGISEI 

will greatly improve our capabilities in 

this area. 
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